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We consider the perturbative (weak field) limit of the small x QCD evolution equation for 
quadrupole, the normalized trace of four Wilson lines in the fundamental representation, which 
appears in di-hadron angular correlation in high energy collisions. We linearize the quadrupole 
evolution equation and then expand the Wilson lines in powers of g Afj, where Afj, is the gauge field. 
The quadratic terms in the expansion (~ g'^ A'^) satisfy the BFKL equation as has been recently 
shown. We then consider the quartic terms (~ A*) in the expansion and show that the linearized 
quadrupole evolution equation, written in terms of color charge density p, reduces to the well-known 
, BJKP equation for the imaginary part of four-reggeized gluon exchange amplitude. We comment 

on the possibility that the BJKP equation for the evolution of a n-reggeized gluon state can be 
^ ' obtained from the JIMWLK evolution equation for the normalized trace of n fundamental Wilson 

Q , lines when non-linear (recombination) terms are neglected. 
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I. INTRODUCTION 



The recent experimental observation of disappearance of the away side peak in di-hadron angular correlation in 
I the forward rapidity region in deuteron-gold collisions at RHIC [l[ has generated a lot of interest in multi-parton 
correlations at high energy (small x). Unlike structure functions in DIS and single inclusive particle production in 
i' hadronic collisions which are sensitive to dipoles (normalized trace of two Wilson lines), di-hadron correlations probe 
1^ , correlators of higher number of Wilson lines {2, 3]. Therefore one has the opportunity, for the first time, to investigate 
these higher correlators experimentally through studies of angular and rapidity correlations in di-hadron production 
cross section in high energy hadronic collisions. Such studies can teach us much about the intrinsic correlations in 
^-H ' the hadronic or nuclear wave functions which are not accessible in single inclusive particle production or in studies of 
^ ] structure function in DIS. 

. Higher correlators of Wilson lines appear in two-hadron production cross section in any dilute-dense collision at high 
' energy where analytic calculations are possible. Classic examples of such asymmetric collisions are proton-nucleus 
collisions (see [J for a review) in the fragmentation region of the proton, and in DIS close to the virtual photon 
remnants ^ . Two-gluon production cross section in DIS has been considered in 0] while two-parton production cross 
section in proton- nucleus collisions has been investigated in 0, d-HI. In all cases, the cross section involves correlators 
of higher (more than two) number of Wilson lines, the most iinportant being the quadrupole operator. Evolution 
equations for these higher point correlators have been derived 0, 0, fiot and approximate analytic expressions for them 
have been developed using a Gaussian model [Tlj and approximate analytic solutions have been proposed [12i |. Very 
recently, powerful lattice gauge theory techniques have been applied to solve the JIMWLK evolution equation which 
then allows a systematic and detailed numerical study of the properties of these higher point correlators • 
5^ ' Here we study the evolution equation for the quadrupole operator in the weak field limit. A first study of this 
has already been performed in (IQ] where it is shown that the quadrupole evolution equation reduces to a sum of 
BFKL equations for the dipole operator in the limit where the dipole is expanded in powers of the gluon field and 
quadratic terms in gluon field are kept. Here, we go beyond the quadratic expansion and show that the quartic terms 
in the expansion of the linearized quadrupole evolution equation satisfy an equation which is identical to the BJKP 
equation [3, [Tsl for the imaginary part of the four-reggeized gluon exchange amplitude. This should be very useful 
since there is an extensive literature on the properties of the BJKP equation which may give us further insight on the 
properties of the JIMWLK equation in the limit where one may ignore non-linear terms. 



^ Particle production in the very forward rapidity region in proton-proton collisions at very high energy falls into this category also. 
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II. QUADRUPOLE EVOLUTION EQUATION 

We start by defining the quadrupole operator Q as 

Q{r,f,s,s) = j^tvVrV^V-sV} (1) 
where Vr = V {rt) is a, Wilson hne in the fundamental representation in the covariant gauge 

and A^{x~ ,rt) = 5{x~)a{rt). The gauge field a{rt) is related to the color charge density via dfa"-{rt) ^ 
gp°'{rt) and r,f,s,s etc. denote two-dimensional coordinates on the transverse plane. The evolution equation for 
the quadrupole was derived in 0| in the large Nc limit and using Feynman diagram techniques. It has been recently 
re-derived using the JIMWLK equation where it was shown that there are no Nc suppressed corrections. Here we 
outline the derivation using the JIMWLK formalism jl6|J where the evolution {y = logl/x) of any operator is given 
by 
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(Pz {x — z) ■ (y ~ z) r J. J. J. 

TT J (27r)2 {x~zf{y~zy + ^^^^ " ' ^^^^ 
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and [/ is a Wilson line in the adjoint representation. The derivation of the quadrupole evolution equation is straight- 
forward but tedious. It involves functional differentiation of the Wilson lines and repeated use of the identity 
[[/(^)]a6^6 = V'<{r)t'' V{r). The result is 
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where the 5" matrix is defined as 



S[r,f) = —ivVrV} 



Q{r,r,s,s) 
S{r,s) S{f, 5) 
^(r,r-) 5(5, s)^J(5) 

(6) 



We will refer to the first four lines in this equation as "real" and the last three terms as "virtual" terms in coordinate 
space. This is to distinguish them from the real and virtual terms in momentum space after we Fourier transform the 
equation since there is no one to one correspondence between the real and virtual terms in coordinate and momentum 
spaces. We have also verified that this equation is exact in the sense that there are no Nc suppressed terms in the 
equation itself (note that models used to evaluate the color averaging denoted by < ■ • • > may introduce sub-leading 
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Nc terms). It also agrees with the previous results for the quadrupole evolution equation 0, [l^j- The S matrix 
satisfies the BK evolution equation [T3| given by 
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(7) 



Unlike the dipole kernel in the BK equation which allows a probabilistic interpretation in coordinate space, the same 
is not true in the quadrupole evolution equation due to terms with negative signs. Even though the individual kernels 
in eq. ([5]) are just the standard dipole kernels [l^, it is still perhaps useful to explain in a more intuitive way, the 
various terms that appear in eq. ([5]). The first four lines in eq. ^ are the "real" corrections and come from the 
third and fourth terms in eq. (|4]). One can rewrite any kernel in eq. ([5]) in a way which may look more familiar and 
facilitates the comparison with the standard dipole emission kernel. For example, the kernel in the first line on the 
right hand side of eq. ([S|) can be written as as 



1 {r — z) ■ {r ~ z) (r — z) ■ {s — z) (r — z) • (s — z) 

(r — z)2 (r — z)2(f — z)2 (r — z)^(s — z)^ (f — z)-^{s — z)^ 



(8) 



with a similar form for all the other kernels. Here the first term corresponds to a gluon being radiated by a quark line 
represented by V(r). If it is absorbed by the same quark line in the amplitude, it leaves the quadrupole unchanged 
and will correspond to a "virtual" correction. On the other hand if it is absorbed by the same quark line but in the 
complex conjugate amplitude (so the gluon line crosses the cut), it will multiply the quadrupole with the coordinate 
r replaced by z and a dipole with coordinates r, z. This will be part of the "real" corrections. The second term above 
corresponds to the case when the quark line, represented by V{r), in the quark anti-quark system represented by 
V{r) and V{f) radiates a gluon with transverse coordinate z. If the radiated gluon does not cross the cut line and is 
absorbed by the anti-quark line at f it becomes part of the "virtual" corrections. On the other hand if the radiated 
gluon at z crosses the cut and is then absorbed by an anti-quark line in the complex conjugate amplitude, it breaks 
the original quadrupole into a quadrupole with coordinate r replaced by z and a dipole at r, z. This is part of the 
"real corrections. All other terms have a similar interpretation. 

To investigate the weak field limit of this evolution and to make our approximations more transparent, it is more 
useful to work with the T matrices, defined as Tq = 1 — Q and T = 1 — S. It is easy to see that all kernels multiplying 
1 (when we switch from Q, S to Tq, T) add up to zero. Therefore, eq. (O is re-written as 
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A. The weak field limits 



It is useful to consider the above equation for Tq in the weak field (dilute) limit where all sizes are much smaller 
than the inverse saturation scale, i.e., \a — b\ << 7^- for any external coordinates a,b. In this limit the non-linear 
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terms {Tq T and T T) in eq. ([9]) may be dropped and we get 
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(10) 

To proceed further, we first consider the two-gluon exchange limit, i.e., the BFKL equation [19]. Since Tq and T 
include multiple gluon exchanges, we need to linearize them, i.e., take the single (reggeized) gluon exchange limit. 
This corresponds to expanding each of the Wilson lines in the definition of Tq and T to first order in the gauge field 
a and then keeping terms of the order a^. In this limit (note the relative sign which appears when taking both a's 
from either V^s or V^^'s rather than taking one a from a V and another a from a V^^) 

TQ(r, f, s, s) T{r, f) + T(s, s) - T{r, s) - T{f, s) + T{r, s) + T{f, s) (11) 

Using eq. ([TT]) in both sides of eq. (fTO|) we get the BFKL equation for each T of a given argument. For example, 
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This limit was already considered in fTo*! and the correspondence with BFKL was shown. We also note that this 
relation still holds when the evolution equation is written in terms of the color charge density p rather than the gauge 
field a. 

The next interesting case is to consider 0{a'^) and see whether our evolution equation reduces to the well-known 
BJKP equation governing the evolution of four reggeized-gluon state in the dilute limit. To do this, again we first 
ignore the non-linear terms in the evolution equation, then we expand the Wilson lines and keep terms of the order 
a* in eq. (|10p. Since the BJKP equation is written in momentum space, we will start by Fourier transforming Tq 
(ignoring T at the moment) to momentum space and disregard any contribution which leads to a vanishing external 
momentum. We define ^ 

Ti{hMMM) = I d^rd^fd^sd^se'^^'-'' + ^^-'' + ^'-' + ^'-'^Ti{r,f,s,s) 



(14) 

where Zi, Z2, ^3, Ia arc two-dimensional external transverse momenta satisfying overall transverse momentum conserva- 
tion so that there are only three independent momenta. This corresponds to having a choice in picking the origin of 
the coordinate space on the transverse plane. One can then right away see that the last term in ([S]) convoluted with 



We will use the notation T4 iiere to denote the ~ 0{a'^) terms in the expansion of Tq so that T4 = Tr [a a a a]. 
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T4{z, f, s, s) will give a S^{li) since it does not depend on coordinate r. A similar argument shows that the last term 
in each kernel in the first 4 lines in eq. (|10l) (the "real" terms) will lead to a delta function which sets one of the 
external momenta to zero. Since the external momenta of the reggeized-gluons are assumed to be finite (non-zero), 
all these terms can be safely ignored. We now consider the contribution of the "virtual" terms, line 5 in eq. ^ 
Upon Fourier transforming, we get 
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with a cyclic permutation of the external momenta in the second term understood. The first term is part of the 
virtual corrections while the second term is part of the real corrections in momentum space. Let us consider now the 
contribution of "real" terms. Fourier transforming the non-zero terms in the first line of eq. (|10l) gives 
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The first term in eq. (|16p is part of the virtual corrections (in momentum space) while the second term is part of the 
real corrections. With a slight rearrangement of the first term one can rewrite the contribution of the first line in eq. 
(HUl) as 
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It is clear that the first term in the square bracket in eq. (IT71) partially cancels the first term in eq. ([T5|) . This 
cancellation becomes complete when we include the similar contributions from the lines 2 — 4 in eq. (ITOl) so that the 
only virtual correction left so far is the second term in the square bracket in (jl7p . Including the contribution of the 
second line to the real part (only the terms which lead to T^^ with the same argument, at the moment) gives 
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where • • • stands for real contributions obtained by appropriate permutation of the external momenta. Finally we 
note that the term proportional to h ■ I2 comes from keeping 0{'^ a^) in the expansion of Vz and setting one of the 
other V^s to unity, for example, taking Vf = 1 and a^iz) in the first line of eq. (llOp. It is clear that the virtual 
terms in eq. ([T8| are already in exact agreement with one gets from B JKP equation |l4l. ITsI l but the real terms look 
different. To show agreement of the real terms with the BJKP equation, we rewrite this equation for color charge 
density p rather than the gauge field a (this does not affect the virtual corrections). To this end, we note that the 
square bracket in the real term in eq. ()18|) can be rewritten as 
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Recalling the relation between gauge field a and color charge density p, 
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and defining 74(^1,12,^3,^4) — Tr p{li)p{l2)p{l3)p{l4), we multiply both sides of eq. with llll^ll which 

effectively removes the external legs. Eq. ((T8)) can then be written as 
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This is our final result and corresponds to the evolution of T4 after one step in rapidity y as depicted (the real part) 
in Fig. (1). We have checked that it agrees with the expressions given in [l5|. 
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h I3 I I2 h 



FIG. 1: Evolution of the four-point function T4, after one step in rapidity as given by eq. (|20|) . Shown is one of the real diagrams 
only and the dashed line represents a cut. 

There are several points that need to be clarified; first, we have completely disregarded the dipole terms T) here 
even though they also contain O(a^) terms. Since T{r,s) depends only on two external transverse coordinates r,s, 
0{a'^) terms will necessarily involve two pairs of gauge fields at the same point. Assuming rotational invariance on 
the transverse coordinate plane, this leads to setting two of the external momenta equal to each other which takes 
one back to the BFKL ladders. Therefore, these terms are not relevant for our purpose. A second point is the color 
averaging denoted by < • • • >. We have not made any assumptions about the color averaging [20| and the evolution 
equation derived is independent of how one performs this averaging. Furthermore, the overall color structure of the 
equation seems to be more general than the BJKP equation since here one has a trace of four color matrices in the 
fundamental representation on both sides of the equation. This trace could be written in terms of products of the 
group structure constants jabc ^^abc -^^j-^ej-eag the BJKP equation is for the exchange of four reggeized-gluon state 
in a symmetric color singlet state. One expects that 5 S terms would lead to a topology which is equivalent to exchange 
of two independent BFKL ponierons which would then be disregarded. Therefore, one would only consider the color 
symmetric structures involving d's. 

In summary, we have shown in this preliminary study that the JIMWLK evolution equation for the quadrupole 
operator can be reduced to the BJKP equation for the real part of the four reggeized-gluon exchange amplitude. 
To do this, we first ignore the non-linear (recombination) terms in quadrupole evolution equation, and then expand 
the Wilson lines in terms of the gauge field (or equivalcntly, the color charge density). This approximation should 
be valid when the external momenta are larger that the saturation scale, i.e., in the dilute region. The quadrupole 
evolution equation reduces to a sum of independent BFKL equations in O(p^) and to the BJKP equation when 
one looks at the terms of order ^ p"*. This suggests that the JIMWLK evolution equation for the n-pole operator 
< TrV{xi) V^X2) ■ ■ • V^{xn) > in the linear limit (dilute region) may be equivalent to the BJKP heirarchy for 
the imaginary part of the n reggeized-gluon exchange amplitude. This would be very useful since there is much that 
is known about the BJKP equation and its properties but not much is known about the properties of the JIMWLK 
equation in analytic form. Proving the equivalence between linearized JIMWLK and BJKP equations may not be so 
difficult since the JIMWLK evolution equation for < TrV{xi) V''{x2) ■ ■ ■ V^ixn) > can almost be written down 
by inspection in analogy with the pattern seen in eq. ([S]). The problem reduces to keeping track of which quark line 
radiates a gluon and counting all the possibilities since all emission kernels are just the standard dipole kernel. It 
would also be interesting to investigate the connection between the non-linear terms in the JIMWLK equation and 
multi-pomeron vertices employed in reggeized-gluon approach to high energy scattering. These issues are beyond the 
scope of this preliminary work and will be reported elsewhere. 
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